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SECTION  1 

INTRODUCTION  AND  STATEMENT  OF  PROBLEM 


In  order  to  optimize  the  effectiveness  of  expensive  weapons,  it  is  essential  for 
intelligence-gathering  systems  to  not  only  detect  and  locate  potential  targets,  hut  also  to 
properly  identify  those  targets  that  are  of  interest.  As  a  result,  the  accurate 
identification  of  military  targets  has  been  a  subject  of  great  interest  in  recent  years. 
The  emphasis  has  mainly  been  on  solving  the  problems  associated  with  separating  target 
returns  from  noise  and  background  clutter  interference,  and  with  the  subsequent 
identification  of  those  targets  that  are  of  interest. 

The  Georgia  Institute  of  Technology,  Engineering  Experiment  Station  (GIT/FES) 
has  been  actively  engaged  in  studies  to  identify  stationary,  military  targets  by 
radar  The  overall  objective  has  been  to  analyze  and  synthesize  radar  signal 

processing  techniques  for  application  aboard  an  airborne  platform  with  stationary, 
tactical  weapons  as  primary  targets.  The  present  report  contains  the  unclassified  results 
of  target  discrimination  efforts  conducted  by  GIT/EES  during  the  period  from  September 
20,  1978  through  December  19,  1981;  the  research  program  is  under  sponsorship  of  the 
U.  S.  Army  Electronics  Research  and  Development  Command  (USA  ERADCOM)  through 
the  U.  S.  Army  Research  Office  under  Contract  Number  DAAG29-78-C-0044. 

In  order  to  better  understand  this  report,  it  would  be  worthwhile  to  briefly 
consider  the  process  of  identifying  a  target.  Referring  to  Figure  1.1,  the  identification 
process  may  be  viewed  as  a  sequence  of  three  steps:  the  detection,  discrimination,  and 
classification  stages.  The  detection  process  is  a  relatively  simple  one  and  is  designed  to 
identify  all  radar  resolution  cells  in  the  search  area  which  contain  potential  targets  while 
keeping  the  number  of  false  alarms  from  clutter  and,  primarily,  from  noise  at  a 
relatively  small  level  (the  probability  of  false  alarm  per  cell  is  typically  10"^  or  less). 
Since  the  number  of  cells  which  must  be  investigated  is  large,  the  detection  procedure 
should  be  relatively  simple  and  fast;  the  techniques  employed  for  this  purpose  are  well 
known  and  will  not  be  discussed  in  this  report.  The  second  step  in  the  identification 
process,  the  discrimination  stage,  provides  a  further  attempt  at  reducing  the  number  of 
false  alarms.  For  this  purpose,  efforts  have  been  directed  toward  the  use  of  techniques 
which  exploit  discriminants  such  as  those  that  may  result  from  frequency  agility,  from 
utilization  of  polarization,  or  both.  The  final  step  in  the  identification  process,  the 
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classification  stage,  consists  of  making  a  decision  with  respect  to  the  identity  of  the 
target  by  assigning  it  to  one  of  various  possible  categories. 

Results  of  a  theoretical  study  to  determine  the  performance  of  one  particular 
class  of  discrimination  algorithms  utilizing  polarimetric  information  is  presented  in  this 
report.  Two  types  of  pseudo-coherent  detection  (PCD)  processors  which  utilize 
polarimetric  phase,  the  relative  phase  between  the  horizontally  and  vertically  polarized 
received  signals,  were*,  analyzed  and  their  performance  determined. 

Also  X-band  and  95  GHz  radar  data  were  collected  on  five  military  targets  under 
this  contract  to  be  used  in  e  computer  simulation  of  polarimetric  detection  processors. 
A  description  of  the  datu  collected  and  the  system  used  to  collect  the  data  is  given  in 
this  report.  Howevei ,  the  interpretation  and  processing  of  this  data  will  be  performed 
under  a  later  contract. 
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SECTION  2 

TARGET-CLUTTER  DISCRIMINATION 


2.1  INTRODUCTION 

A  discrimination  processor  is  the  heart  of  any  radar  system  which  must  detect  and 
classify  targets  in  a  background  of  clutter  and  noise.  The  discrimination  processing  can 
be  more  complex  and  time  consuming  than  that  associated  with  detection  processing, 
since  the  detection  processor  will  have  reduced  the  ceih  requiring  further  investigation 
to  a  small  percentage  of  the  total  number  of  cells  in  the  search  area.  In  general,  the 
false  detection  (false  threat  indication)  probability  does  not  have  to  be  as  small  for 
discrimination  as  it  is  for  detection.  However,  because  of  the  importance  of  detecting 
threat  targets,  the  false  dismissal  probability  should  be  low.  Since  it  is  extremely 
difficult  to  find  a  single  discriminant  which  will  provide  superior  performance  over  a 
wide  variety  of  target  and  environmental  conditions,  it  may  be  necessary  to  combine  the 
outputs  of  several  discriminants  to  obtain  a  satisfactory  target/clutter  decision.  Each 
discriminant  should  take  advantage  of  known  differences  in  scattering  characteristics 
between  targets  and  clutter. 

The  radar  sensor  can  measure  the  amplitude,  time  delay,  frequency  (Doppler)  and 
phase  of  the  scattered  signal  as  functions  of  time.  If  the  operational  situation  allows, 
amplitude,  phase  and  Doppler  may  also  be  measured  as  functions  of  frequency,  aspect 
angle  and  polarization.  Most  discriminants  operate  on  one  or  more  of  the  variables 
mentioned  above.  Considering  the  problem  of  detecting  stationary  targets,  it  is  apparent 
that  Doppler  discriminants  are  not  applicable.  Also,  it  is  assumed  that  discrimination 
will  be  performed  on  the  data  collected  in  the  search  (or  detection)  mode.  These  data 
will  probably  have  low  resolution  with  respect  to  the  target  dimensions  so  that  individual 
target  scattering  centers  will  not  be  resolved.  Thus  time  delay  does  not  appear  to  be  a 
viable  dimension  for  discrimination.  However,  we  are  left  with  the  amplitude  and  phase 
of  the  scattered  radar  signal  as  functions  of  time,  polarization,  aspect  angle  and  possibly 
frequency  as  the  quantities  to  be  used  for  discrimination. 

In  this  study,  combinations  of  amplitude  and  polarimetric  phase,  the  relative  phase 
between  the  horizontally  and  vertically  polarized  received  signals,  were  investigated  for 
discrimination  capability.  Results  of  studies  to  determine  the  performance  of  one 
particular  discrimination  algorithm  called  Pseudo-Coherent  Detection  (PCD)  are 
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presented.  Since  the  PCD-phase  only  processor  is  a  single  range  cell  constant  false 
alarm  process,  a  distinction  between  its  performance  in  a  single  cell  and  in  a  multiple 
range  cell  scenario  must  be  made.  Thus,  a  theoretical  analysis  of  PCD  in  a  single  range 
cell  scenario  (intracell)  is  first  addressed,  and  then  the  theoretical  performance  of  PCD 
in  a  non-uniformly,  range  distributed  clutter  scenario  involving  multiple  cells  is 
determined.  This  latter  type  of  performance  may  be  thought  of  as  intercell 
performance.  In  both  of  these  analyses,  the  performance  of  conventional  square-law, 
single  polarization  channel,  detection  is  used  as  a  reference. 

In  addition,  data  collected  at  Aberdeen  Proving  Ground,  Maryland  under  this 
contract  are  described  in  the  last  section  of  this  chapter.  However,  an  interpretation  of 
these  measured  data  will  not  be  presented  in  this  report. 

2.2  1NTRACELL  THEORETICAL  PERFORMANCE  OF  PSEUDO-COHERENT 
DETECTION 

The  theoretical  performance  of  the  PCD  algorithm  in  a  single  cell  (intracell) 
scenario  will  now  be  addressed.  The  block  diagrams  of  the  two  particular  PCD 
algorithms  under  consideration  in  this  section  are  shown  in  Figures  2.1  and  2.2.  In  Figure 
2.1,  the  IF  signals  in  the  horizontally  (H)  and  vertically  (V)  polarized  receive  channels  are 
mixed  together.  Since  they  are  at  the  same  IF  frequency,  a  video  signal  appears  at  the 
mixer  output  which  can  be  represented  as 

v(t)=  H(t)  ’  V(t)  sin(6(t))  (2-1) 


where 


e(t)  =  <t>H(t)  -  4>v(t)  (2-2) 

The  video  signal  can  thus  be  thought  of  as  the  "dot-product"  of  the  H  and  V  signal 
vectors,  H  (_j  ar>d  V  y.  This  implementation  shall  be  referred  to  hereafter  as 
the  pseudo-coherent  detector/dot  product.  Since  phase  information  is  utilized  in  this 
detector  even  though  the  radar  may  be  non-coherent,  this  processing  technique  is  called 
pseudo-coherent. 
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Figure  2.1  Pseudo-Coherent  Detection-Dot  Product 
algorithm  implementation. 
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Figure  2.2  Pseudo-Coherent  Detection-Phase  Only 
algorithm  implementation. 
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In  Figure  2.2,  a  second  PCD  algorithm  is  shown.  In  this  implementation  the 
horizontally  (H)  and  vertically  (V)  polarized  received  signals  are  each  limited  in 
amplitude  and  are  submitted  to  a  balanced  mixer.  In  effect,  the  limiters  and  mixer 
perform  phase  detection.  However,  the  phase  detected  is  the  relative  phase  between  the 
H  and  V  receive  channels,  not  the  absolute  phase  of  each  channel.  The  mixer  output 
voltage,  which  is  proportional  to  the  sine  of  this  relative  phase,  is  subsequently  submitted 
to  a  threshold  and  a  decision  is  made  relative  tc  a  target  present  or  not  present.  This 
relative  phase  shall  subsequently  be  referred  to  as  polarimetric  phase. 

Equations  describing  the  probability  density  functions  of  the  video  signal,  V(t),  for 
eech  of  these  PCD  algorithms  have  been  derived  from  basic  assumptions  and  are 
presented  here.  The  mathematical  details  are  provided  in  Appendix  A.  The  assumptions 
made  are  as  follows: 

constant  amplitude  (non-fluctuating)  target, 

clutter  signals  in  H  and  V  receive  channels  are  uncorrelated, 

clutter  is  Rayleigh  distributed,  and 

target-plus-clutter  is  Ricean  distributed. 

For  the  dot  product  case,  the  probability  densitv  function  was  not  easily 
derivable.  However,  its  Fourier  Transform,  called  the  characteristic  equation/was 
derived  and  is  as  follows: 


c  (w)  = 


Lxp|^(  j cop  ) /( l+to2/4)J  ’  Exp|"( -£>jo^) /( l+uCM)^ 


(2-3) 


(1  +•  a) 1  /  4 ) 


Where 

p  is  voltage  target-to-clutter  ratio,  >/a?/2  a  and  w  is  dependent  variable  of 
the  equation  (not  frequency). 

In  general,  the  discrete  inverse  Fourier  transform  was  used  to  determine  the 
corresponding  probability  density  function.  This  was  accomplished  on  the  computer 
utilizing  the  fast  Fourier  transform  (FFT)  for  each  different  value  of  target-to-clutter 
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ratio,  o  .  Plots  of  the  clutter-only  and  target-plus-clutter  densities  are  given  in  Figure 
2.3  for  several  different  values  of  p. 

For  some  special  conditions,  the  probability  density  function  is  derivable  in  closed 
form.  For  example,  if  the  target -to -clutter  ratio  is  zero  (clutter-only),  then  the  density 
function  is 


f(v)  = 


-  OQ<  V  <  OO 


(2 -4a) 


n 

where  v  is  the  video  voltage  out  of  the  mixer  in  Figure  2.1  and  a  ^  is  the  clutter 
power.  The  clutter-only  density  function  shown  in  Figure  2.3  was  derived  by  use  of  the 
Fourier  Transform  but  is  also  described  by  Equation  (2-4a).  The  proper  threshold  for  a 
given  probability  of  false  alarm  (PF'A)  can  be  determined  from  Equation  (2-4a)  if  the 
clutter  power  (  O  c  )  is  known  as  follows: 


VT=  -0 


Log  (21 
J  e 


PFA) 


(2 -4b) 


Utilizing  the  probability  density  functions  derived  as  described  above,  probability  of 
detection  as  a  function  of  target-to-clutter  ratio  (TCR)  was  calculated  for  a  fixed 
threshold  level.  The  results  are  presented  in  Figure  2.4  with  the  curve  designated  N=l. 
The  effect  of  pulse-to-puise  integration  was  also  determined  where  the  clutter 
component  for  each  putae  was  assumed  to  be  statistically  independent.  This  was 
accomplished  mathematically  by  convolving  the  proper  probability  density  functions 
numerically  on  the  computer  to  obtain  the  integrated  density  function.  The  resulting 
detection  performance  is  shown  in  Figure  2.4  for  N  =  1  to  32,  where  N  is  the  number  of 
puhes  integrated. 

The  performance  of  the  dot  product  detector  is  also  shown  in  another  way  in 
Figure  2.3.  Here,  the  TCR  required  for  a  90%  detection  probability  and  a  10  ^  false 
alarm  probability  is  3hown  versus  the  number  of  pulses  integrated,  N. 
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gure  2.5  Comparative  detection  performance  of  three  types  of  detectors.  TCR  required  versus 
number  of  statistically  independent  pulses  integrated.  =  0.9,  PrA  =  ,„-4. 


The  detection  performance  of  the  PCD-phase  only  has  also  been  determined. 
Equations  describing  the  probability  density  function  of  the  video  voltage,  V(t),  out  of  the 
phase-detector  shown  in  Figure  2.2  were  derived  for  two  types  of  clutter-only  phase 
distributions.  These  will  be  described  in  the  following  discussion.  It  is  assumed  that  a 
circularly  polarized  wave  is  transmitted;  however,  other  transmit  polarizations  may  also 
be  used  with  this  type  of  detector. 

If  clutter-only  is  present  in  a  range  cell,  and  if  the  clutter  consists  of  many 
approximately  equal-sized  scatterers,  theoretically  the  polarimetric  phase  should  be 
uniformly  distributed  and  sin  4>  would  be  distributed  as  shown  in  Figure  2.6  (a).  On  the 
other  hand,  if  the  object  in  the  range  cell  consists  of  only  one  point-source  reflector  such 
as  a  metal  plate,  corner  reflector,  or  the  like,  the  polarimetric  phase  should  be 
concentrated  near +90°  or  a  sin  <J)  value  of  +_1. 

In  the  data  collected  as  described  in  Section  2.4,  the  polarimetric  phase  does  not 
appear  to  be  uniformly  distributed  for  clutter,  but  appears  to  favor  phase  values  near 
zero  degrees.  Therefore,  a  more  general  mathematical  model  was  derived  to  simulate 
this  observed  tendency  as  well  as  uniformly  distributed  clutter  phase.  It  is  presented  in 
the  following  paragraphs. 

The  following  single  pulse  probability  density  function  for  uniformly  distributed 
polarimetric  phase  was  derived  assuming  a  Rayleigh  distributed  clutter  (amplitude)  model 
and  a  Ricean  distributed  target -plus-clutter  (amplitude)  model.  In  addition,  the  clutter 
in  the  H  and  V  channels  were  assumed  to  be  statistically  independent. 


fj(Z)= 


=  ryr-p"  [ 


1+2  .LA  .  B  .  co8(nco8  XZ) 
n  =  l  n 


Where 


■1  <  Z  <  +1 


r(§  +i)PJ]j/2Exp( -pH) 
r  (n+l) 


M(G  +I,n+1,PH) 
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and  where 


r(-  +l)Py/2  Exp(-pv) 

— 2 - M(-  +1  ,n  +  l  ,py) 

r  ( n  +  1 )  2 


Z  is  the  voltage  out  of  phase  detector, 

denotes  single-pulse  probability  density  function, 

M  (a,  b,  c)  is  the  confluent  hypergeometric  function, 
p  y  is  the  voltage  signal-to-clutter  ratio  in  the  vertically  polarized  receive  channel, 
and 

p  l_j  is  the  voltage  signal-to-clutter  ratio  in  the  horizontally  polarized  receive 
channel. 

When  the  target -to -clutter  ratio  in  both  receive  channels  is  zero,  this  probability- 
density  function  becomes 


(2-6) 

f  (Z)=  —  -1— —  ;  -1  <  Z  <  +1 

and  describes  the  density  function  shown  in  Figure  2.6  (a).  Plots  of  the  density  function 
described  by  Equation  (2-5)  for  several  values  of  target-to-clutter  ratio  is  shown  in 
Figure  2.7.  Density  functions  for  pulse-to-pulse  integration  following  the  PCD-phase 
only  algorithm  were  obtained  by  mathematically  convolving  the  appropriate  density 
function  with  itself.  For  a  large  number  of  pulse  integrations,  the  density  function  tends 
to  be  Gaussian.  Probability  of  detection  was  numerically  calculated  and  is  plotted  versus 
target-to-clutter  ratio  in  Figure  2.8.  The  required  target-to-clutter  ratio  for  a  PD  of 
90%  and  a  PFA  of  10“^  is  plotted  versus  the  number  of  pulses  integrated  in  Figure  2.5. 
An  interpretation  of  these  results  will  be  given  later  in  this  section. 

For  the  non-uniformly  distributed  phase  model,  the  following  probability  density 
function  was  derived  for  the  same  assumptions  as  before. 
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(2-7) 


fia)=F\uH>  •  yvx)dx 

Where 


and  where 

Z  is  the  voltage  out  of  phase  detector, 
fj(Z)  is  the  single -pulse  probability  density  function, 
p  is  the  target -to-ciutter  ratio  either  in  the  H  or  V  channel, 

C  is  a  density  shape  factor,  and 

ERF  ( )  is  the  probability  error  function. 

Equation  (2-7)  represents  the  convolution  of  two  phase  functions  and  provides  the 
density  funclion  for  the  voltage  out  of  a  linear  phase  detector; 

z  -  <|>H  -  4>y  (2-9) 
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The  probability  density  functions  of  <)>  (_j  and  <{>  y  are  given  by  Equation  (2-8).  The 
target -to -clutter  ratio  (  p  )  in  this  equation  is  that  in  the  H  receive  channel  when 
f  ^  ( <J>  f_j)  i®  utilized  in  Equation  (2-7)  and  is  that  in  the  V  receive  channel  when  f  ^  (  <f>  y) 
is  used  in  the  equation.  The  density  shape  factor,  C,  is  used  to  determine  the  shape  of 
the  clutter  density  function  as  illustrated  in  Figure  2.6.  Probability  density  functions 
derived  utilizing  Equations  (2-7)  and  (2-8)  for  several  values  of  target-to-clutter  ratio  are 
shown  in  Figure  2.9  where  a  shape  factor  of  unity  was  selected.  Pulse -to -pulse 
integration  detection  performance  was  obtained  by  convolving  the  respective  density* 
functions  as  described  earlier. 

Detection  performance  for  density  shape  factors  other  than  unity  was  also 
derived.  Interestingly,  the  detection  performance  obtained  for  all  shape  factors  was 
identical  to  that  obtained  with  the  uniformly  distributed  polarimetric  phase.  After  a 
close  examination  of  the  equations  utilized,  it  was  found  that  the  degree  of  shape 
imparted  to  the  clutter  density  function  was  nlso  imparted  to  the  target-plus-clutter 
density  function.  Thus,  for  a  given  threshold  setting  the  relative  areas  under  the  clutter- 
only  and  target-plus-clutter  probability  density  functions  did  not  vary  with  the  shape 
factor.  In  other  words,  the  transformation  being  applied  to  obtain  the  desired  clutter 
density  function  distribution  also  changed  the  shape  of  the  target-to-clutter  density 
function  to  the  same  degree.  Therefore,  the  detection  performance  for  all  modeled 
clutter  distributions  is  given  by  the  PCD-phaae  only  detector  curves  in  Figures  2.5  and 
2.8.  '  > 

Equations  for  the  probability  density  functions  of  the  voltage  out  of  a  square-law 
detector  were  utilized  with  the  computer  to  obtain  the  detection  performance  curves 
shown  in  Figurs  2.11.  Plots  of  the  density  function  obtained  are  shown  in  Figure  2.10  for 
several  values  of  target-to-clutter  ratio.  The  purpose  of  generating  these  curves  was  to 
provide  a  reference  to  which  the  polarimetric  processor  could  be  compared.  Also,  these 
square-law  detector  results  provided  validation  of  our  computer  calculation  procedures. 

A  comparison  of  the  PCD-dot  product,  PCD-phase  only  and  square-law  detection 
performances  is  given  in  Figure  2.5.  For  a  PD  of  90%  and  a  PFA  of  10-<:*y  the  required 
TCR  versus  number  of  pulses  integrated  is  shown.  One  notes  that  the  PCD-dot  product 
processor  requires  from  1.6  dB  to  3  dB  less  TCR  than  the  square-law  detector;  the  result 
of  having  an  additional  channel  of  information.  The  PCD-phase  only  and  the  square-law 
detectors  have  equivalent  detection  performance  for  an  integration  of  60  pulses  or 
more.  However,  the  PCD-phase  only  detector  performance  is  not  as  good  for  less  than 
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Detection  performance  curves  for  square  law 
detector  (ideal  threshold  setting).  N  = 
number  statistically  independent  pulses 
integrated,  PFA  *  10“4. 


Figure  2.11. 


60  pulses  Integrated.  Please  note  that  while  the  PCD-dot  product  and  square-law 
detectors  are  dependent  on  an  Ideally  functioning  threshold  selection  processor  (a 
constant  false  alarm  rate  device)  to  provide  the  performances  shown,  the  PCD-phose 
only  detector  does  not.  Since  the  latter  detector  has  limiters  in  it,  its  clutter  density 
function  is  not  dependent  on  clutter  amplitude  levels.  Thus,  a  fixed  threshold  is  adequate 
to  provide  the  detection  performance  shown  in  Figure  2.5 

The  next  section  of  this  report  addresses  the  question  of  what  detection 
performance  can  be  obtained  for  the  PCD-phase  only  and  square-law  detectors  when  the 
threshold  selection  device  (CFAK)  does  not  function  ideally.  Under  this  type  of 
situation,  one  might  expect  the  PCD-phase  only  detector  to  compare  more  favorably 
with  the  square-law  and  PCD-dot  product  detectors. 

2.3  INTERCELL  PERFORMANCE  OF  PSEUDO-COHERENT  DETECTION 

The  objective  of  the  foregoing  discussion  is  to  provide  a  basis  for  evaluating  the 
intercell  performance  of  the  pseudo-coherent  detection  (PCD)  technique.  In  this  section 
of  the  report,  the  PCD-phase  only  technique  is  compared  to  a  conventional  constant- 
false-  alarm  rate  (CFAR)  detection  algorithm.  Since  the  PCD-phase  only  technique  has 
inherent  CFAR  characteristics,  it  is  desirable  to  compare  it  to  this  type  of  detection 
algorithm.  The  performance  of  the  PCD  and  ampiitude-only  CFAR  algorithms  is 
determined  under  both  a  uniform  clutter  distribution  and  a  non-uniform  clutter 
distribution.  The  particular  non-uniformly  distributed  clutter  scenario  assumed  consists 
of  a  discontinuity  at  a  particular  point  in  range.  On  either  side  of  this  point  in  range  the 
clutter  is  uniformly  distributed,  but  at  different  levels.  The  performance  of  the  PCD  and 
CFAR  algorithms  are  evaluated  for  different  target  positions  relative  to  this  clutter 
discontinuity.  The  worst  case  performance  for  the  amplitude-only  CFAR  algorithm  is 
when  the  target  is  adjacent  to  the  discontinuity. 

2.3.1  THEORETICAL  AMPLITUDE  CFAR  DETECTION  PERFORMANCE 

It  is  well  known  that  the  performance  of  a  fixed  threshold  detector  in  the 
presence  of  an  unknown  clutter  level  is  poor.  Figure  2.12  illustrates  the  probability  of 
false  alarm  as  a  function  of  increase  in  noise  power  density  in  dB  from  design  value.  As 
indicated  in  Figure  2.12  by  the  dotted  lines,  if  the  clutter  level  is  3  dB  larger  than  that 

Q 

assumed,  then  instead  of  obtaining  a  probability  of  false  alarm  of  10"°  a  probability  of 
false  alarm  of  10"^  will  be  achieved.  This  means  that  the  probability  of  false  alarm  is 
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increased  by  10,000  to  1.  Thus  the  need  for  a  method  of  estimating  the  clutter  level  as 
accurately  as  possible.  The  device  used  to  estimate  the  clutter  level  to  achieve  a 
constant  false  alarm  rate  is  called  a  CFAR  algorithm.  Typically,  a  CFAR  device 
averages  the  clutter  found  in  cells  adjacent  to  the  target  cell.  For  example,  a  CFAR 
device,  as  illustrated  in  Figure  2.13,  may  consist  of  10  range  cells,  5  cell3  on  each  side  of 
the  test  cell  under  consideration.  The  clutter  levels  in  each  of  the  10  cells  would  then  be 
summed  and  divided  by  10  to  obtain  an  estimate  of  the  average  clutter  level  surrounding 
the  test  cell.  The  detection  threshold  would  then  be  based  on  this  estimated  clutter 
level.  In  homogeneously  spatially  distributed  clutter,  a  10  cell  clutter  level  estimate 
may  still  b8  in  error  because  only  10  samples  are  used  to  determine  the  average  clutter 
level.  Thus,  a  CFAR  device  with  a  finite  number  of  range  cells  does  not  yield  ideal 
detection  performance.  Figure  2.14  indicates  the  detection  performance  one  would 
obtain  with  a  non -fluctuating  target  in  Rayleigh  clutter  and  where  the  threshed  is  set 
for  a  probability  of  false  alarm  10"^.  For  10  cells  (N  =  10),  one  notes  that  the  detection 
performance  is  about  2  dB  worse  than  that  of  the  ideal  detection  curve.  In  other  words, 
2  dB  more  signal-to-noise  ratio  is  required  with  a  10  cell  CFAR  than  would  be  required  if 
the  clutter  level  were  known.  The  curve  in  Figure  2.15  provides  the  detection  loss  in  dB 
as  a  function  of  the  number  of  resolution  cells  used  in  the  CFAR  device.  Note  that  as 
the  number  of  cells  used  decreases,  the  loss  increases.  For  a  10  cell  device  with  a 
probability  of  false  alarm  10~/4  approximately  2  dB  is  lost  because  of  the  error  in  the 
CFAR  estimate  of  threshold  level. 

When  the  clutter  is  not  homogeneously  distributed  in  range,  but  consist  of  a 
discontinuity  a3  illustrated  in  Figure  2.16,  the  performance  of  the  amplitude  CFAR 
algorithm  deteriorates  noticably.  Figures  2.17(a)  and  (b)  indicates  how  the  detection  and 
false  alarm  proabilities  vary  with  the  magnitude  of  a  clutter  discontinuity  and  also  with 
the  location  of  a  target  relative  to  the  discontinuity  for  a  40  cell  CFAR  algorithm16. 
When  the  target  is  adjacent  to  a  discontinuity  in  clutter  level,  the  CFAR  device 
estimates  the  clutter  level  to  be  half  way  between  the  two  clutter  levels.  Thus,  the 
threshold  is  set  too  high  and  the  probability  of  detection  is  lower  than  desired  while  at 
the  same  time  the  probability  of  false  alarm  is  also  lower.  If  the  target  is  a  few  cells 
from  the  clutter  discontinuity,  the  CFAR  performs  better  because  fewer  of  the  cells  in 
the  CFAR  device  contain  the  higher  level  clutter.  Thus  the  clutter  estimate  and  the 
threshold  setting  is  closer  to  that  desired.  If  in  the  10  cell'example  a  target  is  at  least  5 
range  cells  from  the  clutter  discontinuity,  then  only  one  type  of  clutter  influences  the 
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CFAR  clutter  average  and  threshold  setting.  Thu3,  it  is  only  when  the  target  is  within  5 
cells  or  less  of  the  clutter  discontinuity  (in  the  10  cell  example)  that  the  CFAR  device 
begins  to  perform  poorly.  The  curves  in  Figure  2.18  indicate  typical  detection 
performance  of  a  10  cell  CFAR  device  as  a  function  of  target  position  relative  to  a 
clutter  discontinuity  for  a  fixed  ratio  of  clutter  levels.  As  indicated  before,  as  the 
target  gets  closer  to  the  clutter  discontinuity  the  detection  and  false  alarm  probabilities 
decrease.  However,  if  the  target  is  thought  to  be  in  the  higher  clutter,  the  proper 
threshold  should  be  based  on  the  higher  clutter  level,  but  some  of  the  CFAR  cells  have 
the  lower  clutter  in  their  estimate.  Thus,  in  this  case,  the  threshold  is  set  too  low  and 
the  probability  of  false  alarm  is  much  higher  than  that  desired.  As  the  target  moves 
away  from  the  clutter  discontinuity  so  that  all  of  the  CFAR  range  cells  are  filled  with 
the  higher  level  clutter,  then  the  false  alarm  probability  approaches  the  design  value. 

In  general,  the  amplitude  CF AR  detection  algorithm  will  not  be  able  to  detect 
targets  in  the  clear  near  a  clutter  discontinuity  and  the  algorithm  will  falsely  detect  the 
leading  edge  of  the  higher  clutter  level.  For  example,  the  amplitude  CFAR  algorithm 
does  not  work  very  well  adjacent  to  a  tree  line.  Other  algorithms  which  overcome  this 
weakness  of  the  CFAR  processor  are  desirable.  The  PCD-phase  only  algorithm  seems  to 
provide  good  performance  in  this  type  of  clutter  scenario. 

2.3.2  THEORETICAL  PSEUDO-COHERENT  DETECTION 

The  theoretical  performance  of  the  PCD-phase  only  detector  within  a  single  cell 
(intracell)  has  been  described  and  discussed  in  the  previous  section  of  this  report.  In  this 
section,  the  theoretical  performance  of  the  PCD-phase  only  algorithm  in  a  multiple  cell 
(intercell)  scenario  will  be  addressed.  In  particular,  the  performance  of  PCD  as 
compared  to  amplitude-only  CFAR  will  be  considered.  The  block  diagram  of  the 
particular  PCD  algorithm  under  consideration  in  this  section  is  3hown  in  Figure  2.19.  The 
horizontally  (H)  and  vertically  (V)  polarized  received  signals  are  each  limited  in 
amplitude  and  are  submitted  to  a  balanced  mixer.  In  effect,  the  limiters  and  mixer 
performs  phase  detection.  However,  the  phase  detected  is  the  relative  phase  between 
the  H  and  V  receive  channels,  not  the  absolute  phase  of  each  channel.  The  mixer  output 
voltage,  which  is  proportional  to  this  relative  phase,  is  subsequently  submitted  to  a 
threshold  and  a  decision  is  made  relative  to  a  target  present  or  not  present.  This  relative 
phase  shall  subsequently  be  referred  to  a3  polarimetric  phase. 


29 


Tree  Cutter 


Crass  Cutter 
0 

2 


— i - 1 — i - 1 — r 

-7  -6  -5  -4  -3  -2 


“i - 1 - r — i - r 

-1  +1  +2  +3  * 


-i - r 

•6  +6 


+7 


Range  Relative  To  Clutter  Discontinuity  (Range  Cells) 


Target  Position  (Range) 


Figure  2.18.  Typical  detection  performance  expected  with 
10  cell  CFAR  processor. 
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figure  2.19  Pseudo-coherent  detection  implementation. 
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Figure  2.20  Pseudo-coherent  detection  (PCD) 
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In  the  foregoing  discussion,  it  is  ass.imed  that  a  circularly  polarized  wave  is 
transmitted.  If  clutter  only  is  present  in  ths  range  ceil  under  consideration,  and  if  the 
clutter  consists  of  many  approximately  equal-sized  scatterers,  the  polarimetric  phase  is 
expected  uniformly  distributed  as  illustrated  in  Figure  2.20.  On  the  other  hand,  if  the 
object  in  the  cell  consists  of  only  one  point-source  reflector  such  as  a  metal  plate,  corner 
reflector,  or  the  like,  the  polarimetric  phase  will  be  concentrated  near  one  value  of 
polarimetric  phase.  The  value  of  phase  received  will  depend  on  the  scatterer  type.  Odd 
bounce  scatterers  (e.g.,  trihedrsls)  yield  a  -90°  phase  at  illustrated  in  Figure  2.20  while 
even  bounce  scatterers  (e.g.,  dihedrals)  yield  a  +90°  phase  value.  For  tne  situation 
illustrated  in  Figure  2.20,  a  sing'e  threshold  is  adequate  for  detecting  odd  bounce 
targets.  However,  since  both  even  bounce  and  odd  bounce  targets  are  expected,  a  double 
threshold  as  illustrated  in  Figures  2.21  and  2.22  are  required. 

When  non -uniformly  spatially  distributed  clutter,  as  previously  described,  is 
encountered  by  the  PCD-phase  only  processor,  the  detection  and  false  alarm 
performance  expected  is  as  shown  in  Figure  2.23.  Since  the  polarimetric  phase  threshold 
is  independent  of  clutter  amplitude  levels,  a  change  in  clutter  level  does  not  affect  the 
fa);;e  alarm  probability  (PpA).  Thus  the  PpA  is  the  same  for  each  range  cell  as 
illustrated  in  Figure  2.23  and,  therefore,  the  PCD-phase  only  algorithm  provides  a 
constant  false  alarm  rate.  However,  the  target-to-clutter  ratio  does  change  with  clutter 
level  and,  hence,  the  detection  probability  does  change  across  the  clutter  discontinuity  as 
illustrated. 

One  notes  two  major  advantages  of  PCD-phase  only  over  amplitude-only  CF AR 
detention:  first,  there  is  no  false  detection  of  the  clutter  edge  or  discontinuity;  second, 
tiie  detection  performance  in  the  clear  is  not  affected  by  the  presence  of  the  higher 
clutter  level.  Thus  a  target  in  a  ceil  adjacent  to  a  tree  line  could  be  detected  with  the 
PCD-phase  only  polarimetric  processor  while  it  would  less  likely  be  detected  with  an 
amplitude-only  CF'AR  processor. 

2 A  ABERDEEN  DATA  COLLECTION 

Previous  research  at  GIT/EES  has  shown  that  various  polarimetric  processing 
techniques  can  provide  significant  target-lo-clucter  improvement  for  those  target  aspect 
angles  where  the  horizontally  and  vertically  polarized  components  of  the  reflected  wave 
differ  in  phase  by  about  90  degrees;  however,  there  is  insufficient  data  available  tn 
determine  the  probability  of  a  90  degree  phase  difference  for  a  target  aspect  angle 
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Figure  2.21  Illustration  of  need  for  double  threshold 
in  pseudo-coherent  detection  processing. 
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Figure  2.22  Pseudo-coherent  detection  implementation  with  double  threshold. 
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Figure  2.23  Theoretical  pseudo-cohererst  detection  performance  in  a 
non-uniformly  distributed  clutter  scenario. 
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selected  at  random.  Therefore,  GIT/EES  collected  experimental  data  on  several  tactical 
targets  placed  on  a  turntable  using  GIT/EES  polarization  diversity  radars.  The  major 
purpose  of  the  experiment  was  to  generate  data  at  all  aspect  angles  so  that  amplitude 
and  polarimetric  phase  for  tactical  targets  could  be  characterized  statistically.  It  was 
proposed  that  the  measurements  be  conducted  at  the  U.  S.  Army  Aberdeen  Proving 
Ground  since  targets  of  interest  and  a  suitable  large  turntable  were  available  there.  A 
description  of  the  site,  the  targets  used,  and  the  radar  configurations  employed  in  the 
data  collection  is  given  in  this  section. 

2.4.1  FIELD  SITE 

The  measurements  were  made  at  Field  8,  Aberdeen  Proving  Ground,  Maryland. 
Mr.  Jim  Andrews  was  our  contact  at  Aberdeen,  and  provided  personnel  and  assistance  as 
required  for  the  operation  of  the  target  positioner,  motor-generator  power  supply  and 
cranes.  In  addition,  he  supplied  communications  between  the  equipment  tower  and  the 
target. 

Personnel  from  ERADCOM,  Fort  Monmouth,  New  Jersey  erected  two  closely 
spaced  towers  at  the  field  site.  A  building  approximately  6  feet  by  12  feet  wa3  installed 
on  top  of  the  two  towers  to  house  the  equipment  as  shown  in  Figure  2.24.  The  building 
was  approximately  8  feet  high  and  contained  a  6  foot  by  4  foot  window,  covered  with  a 
plexiglas  sheet.  The  window  was  the  only  access  for  loading  and  unloading  large 
equipment.  Personnel  and  small  equipment  entered  through  a  "trap  door"  located  in  the 
floor  of  the  building.  Steps  were  provided  from  the  ground  to  the  building  floor.  The 
tower  was  70  feet  high  and  was  located  approximately  700  feet  from  the  target 
positioner.  A  view  of  the  target  positioner  from  the  tower  is  shown  in  Figure  2.25  and  a 
close-up  view  of  a  tank  on  the  target  positioner  is  shown  in  Figure  2.26.  All  radar  and 
data  recording  equipment  was  located  on  the  tower.  Data  reduction  equipment  was 
located  in  a  van  at  the  tower  base  for  "quick  look"  data  analysis. 

The  target  positioner  was  a  large,  wood  covered  turntable  capable  of  supporting 
very  heavy  targets.  The  minimum  rotational  speed  of  the  positioner,  to  maintain 
constant  velocity  (prevent  jerking  motion),  is  1  degree  per  second  or  6  minutes  per 
revolution.  Georgia  Tech  installed  an  incremental  encoder  and  a  micro-switch  on  the 
positioner  to  provide  an  angular  position  readout.  The  encoder  was  driven  by  a  lsl  timing 
gear  from  the  main  drive  sprocket.  Output  from  the  encoder  was  a  series  of  2,000  pulses 
per  encoder  revolution  or  greater  than  10,000  pulses  per  target  revolution.  The  micro- 
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Figure  2.25  A  view  of  target  positioner  from  the  radar  measurement  tower. 


i  view  of  a  tank  on  the  target  positioner.  The  plywood 
reduce  radar  reflections  from  positioner  is  shown  in  the 


switch  closed  once  each  target  revolution.  These  two  signals  were  processed  and 
recorded  on  the  data  tape.  A  skirt  of  wood  was  used  to  mask  the  effects  of  the  rotating 
components  of  the  positioner  and  to  reduce  the  magnitude  of  the  radar  return  from  the 
positioner.  A  view  of  the  positioner  and  skirt  is  shown  In  Figure  2.26. 

There  are  two  towers  near  the  location  of  the  target  positioner  and  there  was  an 
interference  problem  associated  with  the  tower  nearest  the  positioner  using  the  150 
nanosecond  pulse.  Therefore,  it  was  necessary  to  reduce  the  transmitted  pulse  length  to 
75  nanoseconds  in  order  to  eliminate  this  source  of  interference.  Also,  trees  located  in 
the  vicinity  of  the  measurement  site  were  removed  to  further  reduce  interference. 

Calculations  based  on  the  known  geometry  of  the  measurement  site  indicated  that 
the  effect  of  multipath  would  not  be  a  serious  problem.  Spacing  between  nulls  generated 
by  the  interference  of  direct  and  indirect  rays  at  9.4  GHz  are  in  the  order  of  15  cm.  This 
close  spacing  compounds  the  problem  of  verifying  that  multipath  effects  are  not  present 
since  a  target  (corner  reflector)  which  is  physically  small  compared  to  15  cm  is 
required.  Such  a  small  target  is  difficult  to  detect  in  the  presence  of  the  target  support 
and  ground  clutter.  However,  an  attempt  was  made  to  define  as  well  as  possible  the 
characteristics  of  the  illuminating  electromagnetic  field  in  the  vicinity  of  the  target 
using  both  large  and  small  targets  (corner  reflectors). 

2.4.2  RADAR  SYSTEMS 

Measurements  were  made  at  both  X-band  and  94  GHz  frequencies.  Several 
transmit/receive  modes  were  used  and  the  data  were  recorded  on  a  FM  magnetic  tape 
recorder.  The  radar  data  were  in  the  form  of  stretched  (box-carred)samples  of  single 
cell  amplitude  and  polarimetric  phase  signals  from  the  targets  of  interest.  Target  aspect 
angle  data  was  simultaneously  recorded  with  the  radar  data. 

2.4.3  X-BAND  RADAR 

Three  configurations  of  the  X-band  radar  were  necessary  to  collect  the  required 
data,  however,  the  transmitter  for  all  three  system  configurations  was  not  changed.  The 
transmitter  is  a  master  oscillator-power  amplifier  (MOPA)  chain  using  a  50  mW  varactor 
tuned  Gunn  oscillator  as  the  basic  source.  The  X-band  signal  generated  in  the  Gunn 
source  is  pulse  modulated  through  two  series  connected  Watkins  Johnson  M14A  balanced 
mix  ere.  The  modulated  signal  from  the  two  balanced  mixers  may  be  varied  from  a 
continuous  CW  signal  to  pulses  as  short  as  15  r»s.  This  signal  is  amplified  in  the  first  TWT 
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to  a  level  of  1  watt.  A  second  pulse  modulated  TWT  amplifier  then  amplifies  the  signal 
to  the  final  value  of  1  kW.  The  Gunn  source  may  be  operated  in  a  fixed  frequency  mode 
or  its  frequency  may  be  changed  on  a  pulse-to-pulse  basis  over  a  frequency  range  of  j-150 
MHz. 

Transmit  modes  were  determined  by  the  arrangement  of  the  waveguide  between 
the  transmitter  and  the  antenna.  Receive  modes  were  determined  by  the  point  at  which 
the  received  signal  is  coupled  from  the  waveguide  into  the  receiver  assembly.  The 
receiver  assembly  is  common  to  all  three  transmitter  configurations  and  is  shown  in 
Figure  2.27.  The  simplest  radar  configuration  is  shown  in  Figure  2.28.  Using  this 
waveguide  arrangement  the  system  is  capable  of  transmitting  signals  with  fixed 
horizontal,  fixed  vertical  or  alternating  (on  a  pulse-to-pulse  basis)  between  horizontal 
and  vertical  polarizations.  The  signals  received  in  the  two  receiver  channels  are  either 
parallel  to  (P  channel)  or  orthogonal  to  (C  channel)  the  transmitted  polarization. 

The  second  X-band  radar  configuration  is  shown  in  Figure  2.29.  This  system  is 
similar  to  the  one  shown  in  Figure  2.28  except  that  a  four  port  90  degree  hybrid,  a  phase 
shifter  and  a  compensating  attenuator  has  been  added.  Power  entering  either  input  port 
of  the  90  degree  hybrid  is  equally  divided  between  the  two  output  ports.  One  output  port 
v/ill  either  lead  or  lag  the  others  by  90  degrees  depending  on  which  input  port  is  activated 
by  the  switchable  circulator.  The  two  signals  are  fed  to  the  orthogonal  ports  of  a  dual 
mode  coupler,  resulting  in  a  transmitted  signal  which  may  be  adjusted  from  a  pure  linear 
polarization  in  a  plane  45  degree  to  the  orientation  of  the  dual  mode  coupler  to  a  pure 
circular  polarization.  Linear  polarization  results  when  the  two  signals  are  in  phase  at  the 
dual  mode  coupler  and  circular  polarization  results  when  the  two  signals  are  separated  by 
90  degrees.  The  receiver  in  this  configuration  is  connected  to  the  waveguide  so  that  the 
parallel  and  cross  channels  represent  the  components  of  the  reflected  wave  which  are 
either  parallel  to  or  orthogonal,  respectively,  to  the  polarization  of  the  transmitted 
wave. 

The  third  configuration  for  the  X-band  radar  is  shown  in  Figure  2.30.  In  this 
sytem  the  waveguide  is  connected  the  same  as  previously  except  that  circulators  have 
been  inserted  at  the  antenna  ports.  This  allows  the  same  polarization  modes  to  be 
transmitted  but  the  receiver  is  capable  of  receiving  only  the  horizontal  and  vertical 
components  of  the  reflected  radar  signal. 
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SWITCHING 

CIRCULATOR 


TO  RECEIVER  BOARD 


Transmit  Modus 

1-  fixed  Horizontal 

2.  Fixed  Vertical 

3.  Alternate  H,  V 


Receive  Modex 

!•  Fixed  Horizontal 

2.  Fixed  Vertical 

3.  Alternate  H,  V 


Figure  2.28  X-band  radar  configuration  one 
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SWITCHING 


PHASE 

SHIFTER 


TO  RECEIVER  BOARD 


Transmit  Modes 

1 .  Fixed  R/L  Circular 

2.  Alternate  R/L  Circular 

3.  Alternate  45  Degree  Linear 

4.  Fixed  45  Degree  Linear 


Receive  Modes 

1 .  Fixed  R/L  Circular 

2.  Alternate  R/L  Circular 

3.  Alternate  45  Degree  Linear 

4.  Fixed  45  Degree  Linear 


Figure  2.29  X-band  radar  configuration  two. 


PHASE 


Transmit  Modes  Receive  Modes 

1.  Fixed  R/L  Circular  1.  Fixed  Horizontal  and  Vertical 

2.  Alternate  R/L  Circular 

3.  Fixed  45  Degree  Linear 

4.  Alternate  45  Degree  Linear 


Figure  2.30  X-band  radar  configuration  three. 
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2.4.4  95  GHz  RADAR 


A  complete  block  diagram  of  the  95  GHz  radar  Is  shown  in  Figure  2.31.  The 
system  is  capable  of  transmitting  either  fixed  or  alternating  vertical,  horizontal,  or  45 
degree  linearly  polarized  signals  as  well  as  fixed  or  alternating  right  and  left  circularly 
polarized  signals.  The  receiver  is  connected  so  that  only  the  horizontal  and  vertical 
components  of  the  reflected  radar  signal  are  detected. 

2.4.5  CALIBRATION  AND  PROCEDURES 

Calibration  of  the  radar  return  was  achieved  by  placing  a  corner  reflector  of 
known  radar  cross  section  in  front  of  the  target  and  recording  the  amplitude  of  the  video 
signal  from  this  range  cell.  The  RF  attenuator  in  the  radar  receiver  was  stepped  in  5  dB 
decrements  over  the  dynamic  range  of  the  signal.  Using  this  technique,  the  calibration 
levels  were  in  units  of  radar  cross  section  and  only  a  range  correction  was  required.  This 
i3  a  reliable  method  of  calibration  provided  the  electromagnetic  field  at  the  target  is  not 
contaminated  with  multipath  or  reflections  from  other  objects  located  in  the  same 
resolution  cell.  To  insure  that  the  field  is  uniform  at  the  calibration  target's  location,  a 
vertical  field  probe  was  used  to  determine  if  multipath  effects  were  present. 

A  separate  calibration  or  reference  signal  is  required  for  each  radar  mode.  Each 
day  a  step  calibration  was  made  to  generate  a  receiver  transfer  characteristic  prior  to 
the  first  target  run.  This  step  calibration  was  made  with  a  diplane  oriented  22.5  degrees 
from  vertical  and  using  both  vertical  and  horizontal  polarizations.  After  the  receiver 
transfer  characteristic  was  recorded  each  data  run  was  preceded  by  a  10  second  sample 
of  the  radar  return  from  an  appropriate  target  or  targets  of  known  radar  cross  section. 
For  this  10  second  sample,  the  radars  were  operated  in  the  applicable 
polarization/frequency /resolution  modes. 

2.4.6  DATA  COLLECTION  PROCEDURES  AND  TEST  MATRIX 

Data  were  collected  for  the  various  radar  modes  shown  in  the  test  matrix  for  five 
military  targets  as  well  as  trees  and  grass.  Targets  which  were  placed  on  the  turntable 
were  rotated  in  azimuth  through  720  degrees.  Each  target  was  placed  on  the  turntable  so 
that  it  was  facing  away  from  the  radar  near  the  zero  degree  position  of  the  turntable.  A 
description  of  four  of  the  five  targets  used  in  this  data  collection  is  given  in  Tables  2.1 
through  2.4.  No  information  was  available  on  the  2  1/2  ton  truck  used  in  the  data 
collection.  The  combinations  of  frequency  mode,  transmit  polarization  mode,  receive 


Figure  2.31  95  GHz  radar  transmitter- receiver  block  diagram. 


TABLE  2-1.  DESCRIPTION  OF  HOWITZER,  LIGHT, 
SELF-PROPELLED,  105MM,  M108 


Weight  Classification 

Weight . 

Length . 

Width . 

Height . 

Fuel . 

Range  .  . . 


Cross  Country  Loaded 
Class  20 

46,  921  pounds 

240-  1/2  inches 

124  inches 

120  inches 
137  gallons 
220  miles 


Speed 


35  miles  per  hour 


TABLE  2.2.  DESCRIPTION  OF  TANK,  COMBAT,  FULL  TRACKED, 
105  MM  GUN,  M60  AND  M60A1 


Weight  Classification . 

Weight  (empty) . . 

Length .  . . 

Width  . . . . . . . 

Fuel . . . 

Engine  (continental): 

Cylinders  . . .  . .  , 

Model  Number . . 

T ype  (both  models) . . 

Transmission,  GMC  (Allison  Division): 
Model  . r , 
Type . 


Cross  Country  Loaded  -  Class  50 


M60 

M60A.1 

95,300  pounds 

97,494  pounds 

320  inches 

325  inches 

144  inches 

128.20  inches 

385  gallons 
(Diesel) 

375  gallons 
(Diesel) 

12 

12 

AVDS  1790-2 

AVDS-1790-2A 

90°,  V-Type,  Air-Cooled,  Diesel 

CD  850-6 

CD-850-6A 

Cross  Drive 

Single  Stage 

Differences  Among  Models: 

The  turrets  are  of  different  design.  Also,  the  external  storage  rack  of  the  M60A1 
extends  completely  around  the  resr  of  the  turret  and  the  lower  portion  is  screened. 


TABLE  2.3.  HOWITZER,  LIGHT,  TOWED,  105MM,  M102 


Classification . . . Standard  A 

Weight  Classification . Empty  -  Class  2 

Componentsi 

Cannon  Carriage  Recoil  Mechanism 

M137MOD  M31  M37 


Length: 

Travel  Position 

Width  . . . 

Weight . 

Cannon: 

Weight . 

Length ...... 

Type  of  Breechblock  .... 

Type  of  Firing  Mechanism . . 

Muzzle  Velocity . . . 

Rate  of  Fire . . 

Estimated  Accuracy  Life  of  Tube 

Carriaqe: 

Height 
Length 

Recoil  Mechanism: 

Type . . . . . M37 

Length  of  Recoil . Variable,  from  30  to 

50  inches 

Equilibrator  . . .  Springs 


47  inches 
176  inches 


900  pounds 
131  inches 
Vertical  Sliding 
Wedge  Manual 
Inertia  Percussion, 
Mechanical 
2,225  FPS 
Max  SPM  -  10, 

Sustained  SPM  -  3 
3,0C0  RD 


262  - 1/2  inch 
76  inches 
3,140  pounds 
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TABLE  2.4.  TRUCK,  UTILITY,  1/4-TON,  4x4,  M151,  M151A1,  M151A2. 


Classification . . .  Standard  8 

M151,  M151A1  .  , . .  Standard  A 

Engine  (Ford,  Army  Design,  4  cyclinders)  . . . .  2805-678-1820 

Transmission  (Selective  Synchromesh) . * .  2520-678-1808 

Weight: 

M151  . . . . . . .  2,273  pounds 

M151A1,  M151A2  . . .  3,600  pounds 

Length . . . . . . .  132.7  inches 

Width . .  . . . .  64  inches 

Height . . .  85  inches 

Tires  . . . .  7.00  x  16 

Battery . .  . . . .  2  (12  V)  24  V  System 

Fuel . . . . .  17  gallons 

Range . . .  300  miles 

Speed . . . . . . .  65  miles  per  hour 


Differences  Among  Models: 

The  M151A1  has  an  improved  suspension  system  over  the  M151.  The  M151A2  has 
a  more  improved  front  and  rear  suspension,  two-speed  electrical  wipers,  manually 
operated  washers,  a  one-piece  windshield,  a  mechanical  fuel  pump,  and  integrated 
exterior  lighting  at  front  and  rear  of  vehicle. 
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polarization  mode,  radar  transmit  frequency  and  range  resolution  utilized  in  this 
measurement  program  are  indicated  in  the  Teat  Matrix  in  Table  2.5. 


TABLE  2.5.  TEST  MATRIX  -  VEHICLE  TARGETS 
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SECTION  3 
SUMMARY 


Results  of  a  study  to  determine  the  performance  of  one  particular  class  of 
discrimination  algorithms  utilizing  polarimetric  phase  have  been  presented.  Since  one  of 
these  processors,  the  PCD-phase  only  processor,  is  a  single  range  cell  constant  false 
alarm  process,  a  distinction  between  its  performance  in  a  single  ceil  and  in  a  multiple 
range  cell  scenario  has  also  been  made.  In  a  theoretical  comparison  of  the  PCD-dot 
product,  PCD-phase  only,  and  square-law  detection  (SLD)  processors,  it  was  assumed  that 
ideal  threshold  settings  were  available  for  the  PCD-dot  product  and  square  law 
detectors.  However,  such  an  assumption  was  not  necessary  for  the  PCD-phase  only 
detector.  Under  this  assumed  situation,  the  PCD-dot  product  detector  requires  from  1.6 
to  3  dB  less  target-to  clutter  ratio  (TCR)  for  the  same  PD  and  PFA  than  the  9quare-law 
detector.  The  3  dB  advantage  exists  for  single  pulse  detection  while  the  1.6  dB 
improvement  is  for  a  very  large  number  of  pulse  integrations.  Also  under  the  above 
assumption,  the  PCD-phase  only  and  square-law  detectors  have  approximately  equivalent 
detection  performance  for  an  integration  of  60  pulses  or  more.  However,  the  PCD-phase 
only  detector  performance  is  not  as  good  for  less  than  60  pulses  Integrated. 

However,  when  the  effect  of  adding  a  CFAR  device  to  a  conventional  square-law 
detector  (SLD)  is  taken  into  account,  the  PCD-phase  only  processor  has  better 
performance  than  the  SLD/CFAR  in  some  clutter  scenarios.  Specifically,  one  notes  two 
major  advantages  of  PCD-phase  only  over  SLD/CFAR  detections  first,  there  is  no  false 
detection  of  a  clutter  edge  or  discontinuity;  second,  the  detection  performance  in  the 
clear  is  not  affected  by  the  presence  of  the  higher  clutter  level.  Thus  a  target  in  a  cell 
adjacent  to  a  tree  line  could  be  detected  with  the  PCD-phase  only  polarimetric  processor 
while  it  would  less  likely  be  detected  with  an  SLD/CFAR  processor. 

X-band  and  94  GHz  radar  data  were  collected  on  five  military  tarets  over  720 
degrees  of  aspect  angle.  This  data  was  collected  to  support  computer  simulation  of 
polarimetric  processors  utilizing  measured  data.  Each  target  was  placed  on  a  turntable 
and  radar  data  was  taken  as  the  target  rotated.  In  addition,  data  was  taken  on  nearby 
grass  and  tree  clutter.  Combinations  of  frequency  modes,  transmit  polarization  modes, 
receive  polarization  modes,  radar  transmit  frequencies  (9.4  and  94  GHz),  and  range 
resolutions  were  utilized.  The  interpretation  and  processing  of  this  data  will  be 
performed  under  a  later  contract. 
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APPENDIX  A 

MATHEMATICAL  MODEL  OF  TWO  DISCRIMINATION  PROCESSORS 


The  mathematical  model  of  a  two  channel  radar  receiver  Is  shown  In  Figure  A-l. 
The  transmitted  signal  is  assumed  to  be 


ST(t  )=Acos  (u>0 1  +  9 ( t)]  ;nT£t  <n(T+i) 


(A-l) 


where 

T 

T 

n 

A 


a) 

o 


e  ( t ) 


is  the  interpulse  period 
is  the  pulse  length 
is  the  nth  pulse  transmitted 
is  the  amplitude  (voltage) 
is  the  transmitted  frequency 

is  the  unknown  phase  associated  with  a  non-coherent  source 


[Sv(t)  '  S  (t) 1 

baseband 

Sin  (4»v(t>  -  <|>H(t)) 

baseband 


Detector  B 
f PCD-Phase  Only) 


Figure  A-l.  Model  of  radar  receiver  and  detection  circuitry. 


The  received  signal  is  assumed  to  he  a  large  sum  of  attenuated  and  time  delayed 
transmitted  signals,  all  from  the  same  range  cell.  This  assumption  is  valid  when  clutter 
is  present  and  consists  of  a  large  number  of  scatterers  in  one  range  cell.  Thus,  the 
received  signal  is  given  by 


(A-2) 


SR(t)  = 


N 

Z 

i  =1 


K . Ac  o  s 
1 


t  +  0(  t )  + 


where 

is  the  attenuation  factor  (a  function  of  range,  reflection  coefficient, 
etc.) 

is  the  number  of  scatterers  in  one  range  cell,  and 

is  the  phase  shift  Introduced  by  each  scatterer  and  is  assumed  to  be 
randomly  distributed  over  0  to  2  it  radians. 

The  received  signal  given  by  Equation  (A-2)  may  be  considered  as  having  in-phase 
and  quadrature  components  as  follows: 


Ki 


N 
+  , 


N 

SR(t)=  Z 
R  i=l 


(KjAC08$.  )C08  jjUot+0(  t^J 


(A-3) 


or  as 


Sr(  t)=X(  t)Co8^U)ot  +  0(  t)J  -  Y(  t )  S  i  n£o)Qt  +  0(  t  )J 


(A -4) 


where 


N 

Y(  t )  =  Z  K.A  Sin0.  (A-5) 

i=l  1  ‘ 
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X(t) 


(A-6) 


N 

=  Z  K.A  cost)) 

i=l  1 

X(t)  and  Y(t)  may  be  considered  as  random  processes  for  both  clutter  only  and  clutter 
plus  target. 

For  the  clutter  only  case,  X(t),  Y(t)  are  each  assumed  to  be  Gaussian  distributed 
with  zero  mean  and  with  variance  (  o  equal  to  clutter  power  (Pc).  They  are  also 
assumed  to  be  statistically  independent.  These  assumptions  are  valid  if  the  range  cell 
under  consideration  contains  a  large  number  of  approximately  equal-sized  scatterers. 

When  a  target  (non-fluctuating)  is  present  with  clutter  in  the  range  cell,  the 
combined  target  plus  clutter  signals,  X(t),  Y(t),  may  be  modeled  as  Gaussian  with  non¬ 
zero  means,  a  and  b,  respectively.  Since  the  amplitude  and  phase  modulation  are  related 
to  X(t)  and  Y(t)  as 


the  in-phase  and  quadrature  means,  a  and  b,  are  related  to  the  means  of  A(t)  and  <J>  (t)  as 
follows: 


(A-10) 


Thus,  specification  of  a  and  b  not  only  defines  the  mean  of  the  amplitude  modulation,  but 
also  defines  the  mean  of  the  phase  modulation. 
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The  received  signal  given  by  Equation  (A-4)  may  represent  the  signal  in  either  the 
horizontal  polarization  channel  or  the  vertical  polarization  channel.  Thus,  the  signals  in 
these  two  channels  may  be  written  as 

(A-ll) 

-  Yy(  t)Sin  Jo)Qt  +  0(  tjj 

(A-12) 

-  Y|_|(  t)sinjcuot  +  0(t)j 


Sy(  t)=Xy(  t )  co s  |a>o t  +  0  ( t)] 
and 

SH(  t )  =Xj_|(  Ocosfu^t  +  et  t)J 


In  Figure  A-l,  two  types  of  detectors  are  presented.  Detector  A  utilizes  all  the 
information  in  each  receive  channel-both  amplitude  and  phase  modulation  information. 
This  detector  is  called  the  Dot-Product  detector  since  it  forms  the  product  of  the  signals 
in  the  two  receive  channels.  Detector  B  in  Figure  A-l  has  amplitude  limiters  in  each 
receive  channe1  preceding  the  mixer.  Thus,  the  output  of  the  mixer  provides  only  the 
cosine  of  the  phase  difference.  Note  that  the  unknown  transmitted  phase,  0  (t),  does  not 
appear  at  the  detector's  output.  Thus,  the  radar  does  not  have  to  be  coherent  when  used 
with  either  of  these  two  detectors. 

In  the  subsequent  text,  mathematical  equations  will  be  derived  to  represent  the 
probability  density  functions  at  the  output  of  both  detectors.  As  will  be  shown,  several 
different  clutter  models  will  be  assumed.  The  justification  for  each  model  will  be  given. 

DETECTOR  A  (PCD  -  DOT  PRODUCT) 

The  output  of  the  mixer  in  detector  A  involves  the  amplitude  of  the  signal  in  each 
channel  as  well  as  the  relative  phase  between  channels,  (  «j»  y  -  ^).  Thus,  the  mixer 

output  is  given  by 

(A-13) 

O(t)  =  rSv/(t)  ’  Su(t)l  .  .  . 

^  V  H  'J  baseband 


where  Sy(t)  and  Sj_j(t)  are  defined  by  Equations  (A-ll)  and  (A-12).  Substituting  thuse 
equations  into  (A-13)  and  simplifying,  one  obtains 
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0(t)=xv(t)  *  ^(O+YyU)  *  yH(t) 


(A-14) 


As  indicated  earlier,  X(t)  and  V(t)  can  be  assumed  to  be  Gaussian  distributed  with 
mean  and  variance  related  to  the  target  and  clutter  levels.  To  obtain  the  probability 
density  function  for  the  mixer  output,  0(t),  the  following  two  probability  relationships 
are  used.  First,  the  product  of  two  Gaussian  random  variables,  Z  =  X*Y,  is  given  by^7^ 


fz(Z)=X 


00 

00 


U)|  *Y  ID 


(A- 15) 


where  fXY  (•»•)  »s  the  joint  density  function  of  the  random  variables  X  and  Y.  Second, 
the  probability  density  function  of  a  sum  of  two  random  variables,  Z  =  X  +  Y,  is 


fZ(Z)=/I  f^CZ-Y^dY 


(A-16) 


where  fXY  (.»•)  is  as  previously  defined.  The  derivation  of  the  density  function  for  0(t)  in 
Equation  A-14  is  given  in  Reference  18  for  a  non-zero  mean  case.  The  expression  for  the 
density  function  of  0(t)  is  extremely  complicated  and  difficult  to  work  with.  However, 
its  Fourier  transform,  called  the  Characteristic  equation,  is  also  derived  and  is  a  much 
simplier  expression.  Thus,  as  described  in  Section  2,  the  Characteristic  equation  for  0(t) 
was  utilized  and  is  given  by^®^ 


C( 


E*p[(jwp)/{l*u,2/4)]  ‘  Exp£(-pV)/(l  +  aj2/4)j  (A-17) 

«)=  - - 


(1  +  0)  /4) 


where  p  is  the  voltage  nignal-to-noise  ratio,  A//?a  Q  8nd  w  is  the  Characteristic 
equation  variable  (not  frequency). 


DETECTOR  B  (PCD  ■  PHASE-ONLY) 

The  output  of  the  mixer  for  this  case  will  only  involve  the  phase  of  the  H  and  V 
receive  signals  since  the  mixer  Is  preceded  by  limiters.  Thus,  the  mixer  output  is  given 
by 

0(  t )  =s  in  J^v(  t )  -$H(  t )}  (A-18) 


As  indicated  in  Equation  (A-8),  the  two  instantaneous  phase  functions,  ^  y(t)  and  ^  |_|(t), 
are  related  to  X(t)  and  Y(t)  as  follows: 


(A-19) 


(A-20) 


where  X(t)  and  Y(t)  are  as  defined  in  Equation  (A-5)  and  (A-6).  The  task  at  hand  in  this 
derivation  is  to  determine  the  probability  density  function  of  0(t)  for  several  clutter 
models. 

UNIFORM  PHASE  CLUTTER  MODEL 

For  the  first  clutter  model,  assume  the  instantaneous  receive  phase  for  each 
channel  is  uniformly  distributed  over  0  to  2  *  radians.  This  model  is  thought  to  be  vaiid 
if  the  clutter  is  composed  of  many  approximately  equal  amplitude  scatterers  in  a 
common  range  cell.  The  derivation  of  the  probability  density  function,  p(*)  ,  follows: 

y  <1 

Let  Z  =  w  »  t^,en  <J>  =  tari  (Z).  If  X  and  Y  are  each  zero  mean  Gaussian  random  variables, 
A  ng) 

then  It  can  be  shown  that'  ' 
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where 


p(Z)  = 


(Z2+C2) 


all  Z 


C  —  o  2  /  o  ^ 

a  i  =  standard  deviation  of  Y  and, 

a  7  =  standard  deviation  of  X. 

The  probability  density  of  $  is  related  to  p(Z)  as*^^ 


p^(4>)  =  sec2<frp2(Z-*>tan$)  ;  1*1  <_  tt/2 

or 

Pd)(4»)=  -  —  ;  |4>|  <tt/2 

Tr(C2  +  t«n2<j>) 

For  the  case  of  a  ^  =  a  2  or  C  =  1,  Equation  (A-23)  becorr.es 

(4> )  =  ~  ;  J<t>|  £*/2 

and  the  phase  is  uniformly  distributed.  For  non-unity  values  of  C,  the  phase 
uniformly  distributed  about  the  zero  value  of  4>  •  The  difference  in  $  v  and 
given  by  the  convolution  of  p  ^  (  4>  y)  and  P  ^  (-  $  14)  as  follows! 


Pa(a)=P(b  (*>V®  P*(-^H} 

where 


(A-21) 


(A-22) 

(A-23) 

(A-24) 
is  non- 

4»  H  i8 

(A-25) 
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In  general,  thio  computation  is  done  digitally  on  a  computer.  However,  when  C  =  1,  it  is 
obvicus  that 


p  (o)=  “  }  Ja|  <*/2  (A-26) 

Since  the  phase  for  this  this  case  is  uniformly  distributed,  the  detector  output  density 
function  is  easily  determined  to  be 


p„<°> 


(A-27) 


where  O  is  the  same  as  CXt),  the  detector  output  voltage.  For  non-unity  values  of 
C,  po(0)  is  determined  with  the  computer  digitally. 

A  non-fluctuating  target  added  to  the  clutter  signal  is  modeled  by  allowing  the 
random  variables  X  and  Y  to  take  on  non-zero  means.  The  density  functions  for  0(t)  then 
becomes  more  complicated.  The  derivation  follows  with: 


Z  = 


Y 

X 


PX(X)=  --Li~  Exp  Q-(X-a)  2/2CTl2]J;  all  X 
V 2iro  ^  2 


and, 


Py(y)=  Fxp  f-(Y-b)2/2o22];  al 


/ 2iro  2 


1  Y 


then 


(17) 


(A-28) 


(A-29) 


PZ(Z)=2  J*  Ypx(ZY)py(Y)dY 
00 


(A-30) 
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Substituting  Equation  (A-28)  and  (A-29)  into  (A-3G)  and  completing  the  square  in 
the  exponent  yields 


Pz(Z)=Exp  4  — - 


J  f ,. 

2(o22Z  2+o12)  /  tto^o2 

"  •*  -oo 


(AV-B) 


\  2 


2  2 
2ol  °2 


(A-31) 


dY 


where 


A  =  v^22zZ+ct12  ^ 


B  = 


2  2 
o2  aZ+Oj  b 

/  2~2~  P 

>/ a2  Z  +0^ 


(A-32) 


(A-33) 


The  value  of  the  integral  in  Equation  (A-31)  is  determined  by  using  the  identity  in 
Reference  (19)  with  the  proper  substitutions.  Thus, 


P7(Z) =Exp 


r*  i 
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CM 
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N 
-O 
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1  al°2 

I  2(a  2z2+a2) 

4 

jir(a72Z2+o1 2, 

L  L  1  J 

L 

Exp 


( o22aZ+a*  b2) 


„  2  1,  2 -.2  2v 

Zo  ^o  2(a  2Z  +a  j)  j 


(A-32) 


(a2  aZ+a^  b) 

f  o9  aZ+a. d  \ 

1+ERF  I  1  _  | 

p7(o22z2.ai2)’'2 

V fV2 

where 


a  is  the  mean  signal  level  of  the  in-phase  component  (with  non-fluctuating  target 
present), 


b  is  the  mean  cignai  level  of  the  quadrature  component  (with  non-fluctuating  target 
present),  and 


2  2 

°  1  *  a  2  01,0  c^utter  tev®*8  1°  fc^e  two  orthogonal  components  of  the  signal. 


It  is  convenient  at  this  point  in  the  derivation  to  redefine  some  of  the  parameters 
in  Equation  (A-32).  First  define  a  voltage  target-to-clutter  ratio  for  both  the  in-phase 
and  quadrature  components  as  follows: 


6  ,  -  ■/  R5! 

8q  =  b/  J3>2 


(A-33) 


Then,  define  a  clutter  distribution  shape  factor  as 


C  = 

a 

2 


(A-34) 


As  is  explained  in  Section  2,  when  the  clutter  power  or  variance  in  the  in-phase  and 
quadrature  components  of  the  receive  signal  are  unequal,  then  the  resultant  phase 
probability  distribution  is  not  uniform  but  shaped.  With  the  newly  defined  variables,  the 
expression  in  Equation  (A-32)  becomes 


PZ(Z)= 


C  Exp[-(8I2+8Q2)]  C(8jZ+cBq) 


(Z2+C2) 


FuKc2)512 


(A-35) 


*  Exp 


(6qZ-8ic) 


(Z2+C2 


C)2  1  /eiZ+C0Q  \ 

)  J  L  '  Jc2+Z2  t. 


The  instantaneous  phase  in  either  of  the  two  receiver  channels  (H  or  V)  is  given  by 


♦  =  tan"*(Z) 


(A-36) 
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The  probability  density  for  ^  is  related  to  that  of  Z  as  follows: 


f  ($)=  ( sec  $  )  pz(Z.tan*) 


Making  this  transformation,  Equation  (A-35)  becomes 


(A-37) 


f  (J'-Caec^  ) 


C  Exp  £~(pj  +0q  C(p  j  t an^ +C^q) 

(C2  +  tan2$)  JrrVc^  +  tan^ ) 3^2 


lQtnn^  -f  ,C)21  _  r  /  Bjtan**C69\  _ 

(C2  +  tan2<j>)  J  £  \^C2  +  t  an2$  /_ 

j 


(A-38) 


If  the  polarimetric  phase  of  the  target  return  assumes  a  particular  value,  say  90°,  then 
Equation  (A-38)  is  simplified  by  allowing  the  in-phase  target-to-clutter  ratio,  0  j,  to 
have  zero  value.  Equation  (A-38)  then  becomes 


]  C  Exp  (6‘)  c  B  Exp.  , 

f  (<j>)  =  (3ec  )  < - + - *•- - 

I  (C2  +  tan2(|))  p(C2  +  tan2<j>)3/2 


*6.  t  fin  2^__  | 

.  .C^  +  t  an  2<{>J 


I  1+ERF  # - 

L  Ijsrr 


(A-39) 


To  propertly  normalize  the  above  equation,  one  must  divide  it  by  +  ERF  (  0  )|.  This  is 
the  expression  presented  in  Section  2  as  Equation  (2-8). 

For  the  special  case  where  the  shape  factor  is  unity,  the  probability  density 
function  of  the  detector  output  given  by  Equation  (A-27)  can  be  derived  without  using  the 
computer  to  provide  digitlal  convolution.  The  resulting  density  function  is  presented  in 
Section  2  as  Equation  (2-6). 
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